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基于 mtDNA 和 rDNA 基因 序列 的 中 国 按 蚊 属 
塞 蚊 亚 属 种 类 的 系统 发 育 研 究 
(MBA: HF) 


(第 二 军医 大 学 病原 生物 学 教研 室 ， 上海 200433 ) 


摘要 :【 目 的] 应 用 mtDNA 和 rDNA 基因 特征 重建 中 国 按 败 属 塞 蚊 亚 属 已 知 种 类 的 系统 发 育 关系 , ARER 
种 的 亲缘 关系 。【 方 法 】 对 采 自 中 国 的 按 败 属 塞 蚊 亚 属 Anopheles (Cellia) 20 种 蚊 的 mtDNA-CO II fll rDNA-28S-D3 Jy 
列 进行 测定 和 分 析 ， 以 按 蚊 属 按 蚁 亚 属 Anopheles ( Anopheles ) 的 中 华 按 蚊 An. (An. ) sinensis 和 赫 坎 按 蚊 An. (An. ) 
hyrcanus 为 外 群 , 采用 CO IM D3 单 基 因 , DL" CO TI + D3” 联 合 数据 组 以 邻接 法 (NJ) 、 最 大 简约 法 (MP) 、 最 大 似 
然 法 (ML) 和 贝 叶 斯 法 (BI) 等 重建 这 些 种 类 的 系统 发 育 树 。 [2528] mtDNA-CO IL F rDNA-28S-D3 序列 的 长 度 范围 分 
别 为 685 bp 和 375 ~410 bp, 在 塞 蚊 亚 属 蚊 种 间 的 遗传 距离 分 别 为 0.015 ~0.117 和 0.003 ~0.111。 各 系统 树 显 示 外 
群 被 合理 分 开 , 除 在 COI 树 中 新 塞 蚊 系 为 并 系 外 , 各 系 均 聚 为 单 系 群 , 新 迈 蚊 系 和 迈 蚊 系 亲缘 关系 最 近 。 联 合 数据 
组 构建 的 系统 合意 树 显 示 中 国 塞 蚊 亚 属 各 蚊 种 形成 4 x , ER TN R a iy SEM 蚁 种 团 未 聚 为 单 系 群 外 ,其 他 各 种 
团 和 复合 体 成 员 种 均 分 别 聚 在 一 起 , 各 分 支 的 置信 值 均 大 于 50% 。 【结论 】 本 研究 获得 的 分 子 系统 发 育 树 清 楚 地 显 
示 了 中 国 按 败 属 塞 蚊 亚 属 各 种 类 及 系 之 间 的 系统 发 育 关系 ， 对 其 分 类 和 | 防治 研究 具有 参考 价值 。 
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Phylogenetic relationship among the species of Anopheles subgenus Cellia 
(Diptera: Culicidae ) in China: inferred by mitochondrial and ribosomal 


DNA sequences 

WU Jing, MA Ya-Jun* , MA Ying ( Department of Pathogen Biology, Second Military Medical University , 
Shanghai 200433 , China) 

Abstract. [ Aim] To reconstruct the phylogeny of the known species of Anopheles subgenus Cellia from 
China based on mtDNA-CO 由 and rDNA-288-D3 sequences, and to interpret the interspecific genetic 
relationship among them. [ Methods] mtDNA-CO I and rDNA-28S-D3 genes of 20 species in subgenus 
Cellia and 2 species in subgenus Anopheles of genus Anopheles were sequenced and analyzed. The trees of 
neighbor joining ( NJ) , maximum parsimony (MP), maximum likelihood (ML) and Bayesian inference 
(BI) of these species were constructed based on CO II, D3 and CO [J with D3 conjunction dataset, with 
Anopheles (An. ) sinensis and An. (An. ) hyrcanus as the outgroup. [ Results] The length range was 685 
bp for mtDNA-CO If and 375 -410 bp for rDNA-28S-D3, and the level of pairwise genetic distances among 
Cellia mosquitoes was 0.015 — 0. 117 for CO || and 0.003 -0. 111 for D3 gene. All the phylogeny trees 
showed that the outgroup was divergent. The Neomyzomyia series was monophyletic , but paraphyletic in CO 
Il tree. The relationship between Myzomyia and Neomyzomyia was close. The 50% consensus tree by CO 
H with D3 conjunction dataset showed that the 20 species in subgenus Cellia from China were clustered to 
4 major clades, and all members of each series and group were respectively clustered together with above 
50% bootstrap value, except Anopheles ( Cellia) pseudowillmori not clustered within the maculatus group. 
[ Conclusion] This molecular phylogeny clearly showed the interspecies and inter-series genetic relationship 
among Chinese subgenus Cellia mosquitoes, and could be useful for the taxonomy and control measures. 
Key words: Anopheles; Cellia; molecular phylogenetics; mtDNA; rDNA; phylogeny 
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塞 蚊 亚 属 Cellia s&J8 T GH H (Diptera), WP} 
( Culicidae ) ， 按 蚊 亚 科 ( Anophelinae ), $E "hx JE 
Anopheles。 中 国 已 记录 的 塞 蚊 亚 属 近 30 种 ( 陆 宝 
Bt, 1997; Harbach and Howard, 2007) , 其 中 部 分 种 
类 可 以 传播 严重 危害 人 类 健康 的 疾病 。 系 统 发 育 研 
究 是 从 蚊虫 种 群 的 演化 历史 和 种 系 发 生 规律 分 析 环 
境 变迁 因子 对 其 生物 学 特性 的 制约 和 有 影响 ,是 阐明 
或 解释 现存 蚊 类 的 生态 特性 和 传 病 关系 的 重要 基 
础 。 近 年 来 分 子 标志 的 应 用 推动 了 该 领域 研究 的 进 
一 步 发 展 ( Norris, 2002; Harbach , 2007) , 在 依据 核 
糖 体 DNA (Sharpe et al., 2000; Garros et al., 2005; 
Torres et al., 2006; Hasan et al., 2008) 或 线粒体 
DNA ( Krzywinski et al., 2006; Dusfour et al., 2007; 
Paredes-Esquivel et al., 2009) 的 基因 序列 重建 按 蚊 
属 多 个 种 团 或 复合 体系 统 发 育 关 系 的 基础 上 , ME 
蚊 亚 属 的 系 和 部 分 种 类 的 分 子 系统 发 育 关 系 已 进行 
了 探讨 (Sallum et al., 2002; Manguin et al., 2008; 
Mohanty et al., 2009; Morgan et al., 2009), SHEA 
ees XT HH Es] d EX 2k WC Anopheles hyrcanus FH HI ( Ma 
and Xu, 2005 ) 和 多 斑 按 An. maculatus $8 B] ( Ma 
et al., 2006) 成 员 种 的 分 子 系统 发 育 关 系 进行 了 重 
建 。 为 阐明 中 国 按 蚊 种 间 的 亲缘 关系 , 并 为 蚊 媒 防 
制 策 略 制定 提供 重要 理论 依据 , ASOT SE OT TUS 28 
蚊 亚 属 部 分 种 类 的 核糖 体 DNA 28S 编码 区 第 3 255 
构 域 (ribosomal DNA 28S coding region third domain, 
rDNA-28S-D3) ) 和 线粒体 DNA 细胞 色素 氧化 酶 亚 单 
f II (mitochondrial DNA cytochrome oxidase subunit 
Il, mtDNA-CO IL ) 基 因 序 列 进 行 了 测定 和 分 析 , 并 
重建 其 种 间 的 系统 发 育 关 系 , 探讨 各 分 子 标志 的 信 


息 贡 献 量 。 
1 材料 与 方法 


1.1 材料 

研究 所 用 蚊虫 的 采集 方法 为 诱 蚊 灯 或 人 工 捕捉 
成 蚊 ， 处 死 后 目 然 干燥 ， 单 头 置 离心 管 融 回 实验 
室 ; 种 类 的 鉴别 依据 陆 宝 肇 等 (1997 ) 。 本 研究 涉及 
种 类 和 标本 信息 如 表 1。 
1.2 目的 基因 的 扩 增 和 测序 

T£ Collins 等 (1987 ) 的 方法 , He HZ p. We FE Al ZH 
DNA, # 5$ mtDNA-COI 和 rDNA-28S-D3 基因 的 引物 
参照 文献 设计 、 合 成 (Beard et al., 1993; Sharpe et al., 
2000) (422), 反应 体系 相同 , 均 含 2 pL DNA 模板 、 
1 x PCR 缓冲 液 、0.1 mmol/L dNTPs, 1.7 mmol/L 
MgCl,、2.0 U DNA Tag 聚合 酶 和 0.1 pmol/L 成 对 引 


物 。 反 应 条 件 为 94% 2 min 后 ,以 94%C 30 s. 45 
1 min, 72C 1 min 共 39 个 循环 , 72°C 延伸 8 min, 17 
增产 物 经 1% BEE A 5%o Coldview) 凝 胶 电泳 检测 
Ja, 割 胶 回 收 PCR 产物 , 四 色 获 光标 记 的 双 脱 氧 链 
终止 法 测序 (ABI 3730 测序 仪 ,上海 英 骏 生物 技术 有 
限 公司 ) , 测序 引物 同 扩 增 引物 , 采用 正 、 反 向 测序 
以 保证 序列 的 完整 性 。 

表 2 扩 增 mtDNA-COI 和 TDNA-28S-D3 的 引物 序列 
Table 2 Primers for mtDNA-CO H and 
rDNA-28S-D3 amplification 
目的 基因 Gene 
mtDNA-CO II 


8 [JJ FEF (5'—3^) Primer sequences 
Forward: TCTAATATGGCAGATTAGTGCA 
Reverse: ACTTGCTTICAGTCATCTAATG 
rDNA-28S-D3 Forward: GACCCGTCTIGAAACACGGA 


Reverse: TCGGAAGGAACCAGCTACTA 


1.3 数据 统计 与 分 析 

测序 结果 以 Bioedit 5.0.9 软件 包 进行 核对 ， 必 
要 时 手工 调整 ,确保 序列 的 准确 性 。 采 用 Clustal X 
1.8 比 对 所 获 序列 ,参数 为 默认 值 , 将 比 对 结果 输 
和 人 MEGA 3.1 软件 包 进 行 序列 组 成 分 析 , 包括 变异 
位 点 (variable sites)、 人 简约 信息 位 点 (parsimony 
sites ) Al Dk FE ZH JV ( nucleotide 
composition) 等 ,依据 无 脊椎 动物 的 线粒体 密码 子 在 
MEGA 中 推测 氨基 酸 序 列 , 并 计算 各 蚊 种 间 的 未 矫 
IE P 遗传 距离 。 

在 PAUP 4.0b10 和 MrBayes V3. 1 软件 包 中 进 
行 系统 发 育 分 析 。 分 别 采 用 邻接 法 (neighbor 
joining，NJ) 、 最 大 简约 法 (maximum parsimony, 
MP) 和 最 大 似 然 法 (maximum likelihood, ML) 在 
PAUP 中 构建 系统 树 。 构 建 ML 树 之 前 ， 应 用 
Modeltest 3. 06 软件 包 选 择 其 最 优 进化 模型 。 应 用 
MrBayes 构建 贝 叶 斯 (Bayesian inference，BI) 树 ， 以 
随机 树 为 起 始 树 ,同时 运行 4 个 马尔 可 夫 链 ( Markov 
Chain Monte Carlo, MCMC) , 43247150 000 代 , 每 10 
代 抽 样 1 次 ,用 似 然 值 达到 饱和 之 后 的 数据 获得 贝 
叶 斯 树 。 对 上 述 4 种 方法 构建 的 系统 树 在 PAUP 中 
进行 SH 检验 (Shimodaira-Hasegawa test) , 检验 各 树 
的 拓扑 结构 差异 是 否 显 著 , 对 差异 不 显著 的 系统 树 
构建 50% 合 意 树 。 

本 人 研究 进行 系统 发 育 分 析 时 分 别 依据 各 单 基因 
数据 和 基因 联合 数据 组 , 在 分 析 前 对 联合 数据 组 进 
fT ILD 检验 (incongruence length difference test) ， 检 
验 联合 分 析 的 数据 是 否 具 有 相合 性 。 
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2.1 基因 特征 

本 人 研究 获得 22 种 按 蚊 的 mtDNA-COIT 和 TDNA- 
28S-D3 序列 , 包括 按 蚊 属 赛 蚊 亚 属 的 20 FH, 以 及 
P MEEHRCBOI Je B PEPE ly A PR iy. PA PSE A RF 
征 如 下 。 
2.1.1 mtDNA-CO II : 对 本 研究 中 22 种 按 蚊 的 
mtDNA-CO II 序列 进 行 比 对 ,结果 显示 mtDNA-CO 
ak A 685 bp， 为 一 个 完整 的 阅读 框 ， 编 码 由 
228 个 氨基 酸 残 基 组 成 的 CO H ,其 起 始 密码 子 为 
ATG, 以 一 个 了 作为 潜在 的 终止 密码 子 , 未 发 现 碱 
基 的 插入 和 缺失 。 

在 685 个 比 对 位 点 中 内 群 20 种 存在 191 个 变 
异 位 点 ,150 个 简约 信息 位 点 ，A + 了 T 平 均 含量 为 
74.4% 。 在 密码 子 的 不 同位 置 ， 碱 基 组 成 差异 较 
K, 其 中 AT 含量 在 第 3 位 为 92.7% ,高 于 第 1 位 
的 63.6% 和 第 2 位 的 66.8% 。228 个 氨基 酸 残 基 比 
对 的 结果 显示 , 保守 位 点 214 个 , 占 93.86% ; 变异 
位 点 14 个 , 46.14%. 

各 蚊 种 间 的 遗传 趾 离 显示 ,内 群 与 外 群 间 范围 
为 0.088 ~0.115, 内 群 蚊 种 间 为 0.015 ~0.117, 其 
中 大 劣 按 蚊 与 大 劣 按 蚊 D 之 间 的 核 苷 酸 变异 最 小 ， 
大 劣 按 蚊 与 塞 沃 按 蚊 间 变异 最 大 ( 表 3)。 
2.1.2 rDNA-28S-D3: 本 研究 22 种 按 蚊 的 rDNA- 
28S-D3 序列 长 度 范 围 为 375 ~ 410 bp, 校对 排 齐 后 


显示 ,内 群 中 核 苷 酸 对 比 位 点 共 446 个 , 其 中 变异 
位 点 165 +, 简约 信息 位 点 135 +, T, C, AMG 
的 平均 含量 为 18.7% , 27.7% , 23.2% Fil 30.4% 。 

各 米 种 间 遗 传 距离 的 计算 结果 显示 ， 内 和 群 与 外 
群 间 范 围 为 0. 108 ~0.149, 内 群 种 类 间 的 核 车 酸 变 
异 范 围 为 0.003 ~0.111, 其 中 大 劣 按 星 与 大 劣 按 蚊 
D、 多 王 按 蚊 与 达 罗 上 毗 按 蚊 间 序列 变异 最 小 ; 腹 簇 
按 娄 与 伪 威 氏 按 蚊 间 变 异 最 大 ( 表 3) 。 
2.1.3 “COI + D3” 联 合 数据 组 

本 研究 对 获得 的 CO ILRI D3 数据 进行 组 合 , 连 
接 后 得 到 “COI + D3 ”联合 数据 组 ,该 组 的 核 苷 酸 
比 对 位 点 共 1 129 个 , 其 中 变异 位 点 369 个 , 简约 
信息 位 点 298 个 。 对 “CO II + D3” 联 合 数据 组 的 数 
据 进 行 ILD 检验 , 结果 显示 两 组 数据 具有 不 相合 性 
(P 20.001), 
2.2 系统 发 育 树 

本 研究 分 别 依 据 CO I, D3 和 “CO + D3” X 
合 数据 组 , 采用 MP 法 、NJ 法 、ML 法 和 BI 法 构建 
系统 发 育 树 。 应 用 Modeltest 对 上 述 3 组 数据 的 序 
列 比 对 结果 分 别 进行 等 级 似 然 比 检验 ,结果 显示 最 
适合 的 DNA 进化 模型 均 为 CTR +I + G 模型 。 
2.2.1 mtDNA-CO II : 以 中 华 按 蚊 和 赫 坎 按 蚊 为 外 
BÉ, 依据 COI 单 基因 序列 构建 了 4 种 系统 发 育 树 , 
SH 检验 的 结果 显示 ML 树 似 然 值 最 小 , 且 与 MP 
BJ. NJ 树 和 BI 树 的 似 然 值 差 异 显著 (P <0.05), 
提示 ML 树 能 较 好 解释 CO 了 数据 ,显示 外 群 被 合 
理 分 开 ， 内 群 形成 5 小 支 ( 图 1)。 
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图 1 依据 mtDNA-COI 序 列 构建 的 塞 蚊 亚 属 ML 树 (CTR +7T+C 模 型 ) 


Fig. 1 


The ML tree of Subgenus Cellia mosquitoes by mtDNA-CO JI dataset (GTR + I + G model) 
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2.2.2 rDNA-28S-D3: 本 研究 对 依据 D3 单 基 因数 
据 构 建 的 4 种 系统 发 育 树 进 行 SH 检验 , 结果 显示 
BI 树 似 然 值 最 小 , 与 MP 树 的 似 然 值 差 异 显著 (P < 
0.05), 5 NJ WA ML 树 的 差异 不 显 车 (P >0.05)， 
提示 NJ, ML 和 BI 树 都 能 较 好 解释 D3 数据 ,保留 


100 
1.00 


58 
0.85 


三 者 共同 的 分 支 , 建立 50% 合意 树 ( 图 2) 5 ERR 
显示 外 群 种 类 与 内 群 分 开 , 内 和 群 种 类 形成 4 文 , NJ 
树 各 分 文 的 置信 值 大 于 58% , BI 树 的 后 验 概率 在 
0.851 Eo 


94 maculatus 
1.00 dravidicus 
willmori 
sawadwongporni 


seudowillmori . . 
p Neocellia serles 


pattoni 
stephensi 
splendidus 
philippinensis 
nivipes 
annularis 
99 minimus 
Bho harrisoni 





aconitus Myzomyia series 
culicifacies 

100 dirus 
dirus D Neomyzomyia series 
tessellatus 





vagus ee Pyretophorus series 
sinensis 
hyrcanus 





Outgroup 


图 2 依据 rDNA-28S-D3 序列 构建 的 塞 蚊 亚 属 NMJ,， ML 与 BI 的 50% 6 ER 
Fig. 2 The 50% consensus tree (NJ, ML and BI) of Subgenus Cellia mosquitoes by rDNA-28S-D3 dataset 
分 支 之 上 数据 为 NJ 树 的 置信 值 , 之 下 为 BI 树 的 后 验 概 率 。NJ bootstrap and BI posterior probability values are shown above 


and below each node ，respectively. 


2.2.3 “COI +D3” 联 合 数据 组 : 对 依据 “COI + 
D3” 数 据 组 构建 的 塞 蚊 亚 属 ML 树 似 然 值 最 小 , 与 
BI 树 的 似 然 值 差 异 显著 (P<0.05), 与 NJ 树 和 MP 
树 差异 不 显 车 ( 忆 >0.05) , 提示 NJ, MP 和 ML 树 都 
能 较 好 解释 联合 数据 组 ， 故 建立 50% 合意 树 
(图 3) 。 合 意 树 显示 外 群 被 合理 分 开 , 内 群 种 类 形 
M4 x: 由 新 塞 蚊 系 11 种 : (( 威 氏 按 蚊 ， 塞 沃 按 
B. REOR, AR ME) OBITU R 
B) SC, JEEP) 环 纹 按 蚊 ,斯 氏 按 蚊 ， 
美 彩 按 蚊 ) ; @ 新 迈 蚊 系 4 种 : (CK RC, KE 
按 蚊 D) EFL Oe, 棋 斑 按 蚊 ); QA 4 种 : 
((( 微 小 按 蚊 , 哈里 森 按 蚊 ) 马 涉 按 蚊 ) 库 态 按 蚊 ); 
@ 带 热 蚊 系 : PEE DC, NJ 树 和 MP 树 的 置信 值 显 
示 外 和 群 与 内 群 分 离 支 持 度 高 ， 内 群 中 分 支 必 和 中 置 
信和 值 均 大 于 96% , 但 分 支 由 中 有 的 小 于 30% 。 


3 讨论 


本 研究 重建 中 国 塞 蚁 亚 属 蚊 种 的 系统 发 育 关系 


依据 的 分 子 特征 包括 线粒体 和 核糖 体 DNA 的 基因 
序列 , 虽然 CO II A D3 基因 序列 突变 频率 不 同 , D3 
比 CO IL SF, 但 在 本 研究 的 蚊 种 中 均 显 示 出 良好 
的 种 特异 性 , 与 文献 报告 的 在 其 他 蚊 种 中 的 应 用 一 
$i ( Sallum et al., 2002; Mohanty et al., 2009; 
Morgan et al., 2009) ,应 是 重建 塞 蚊 亚 属 蚁 系统 发 
育 关 系 理想 的 分 子 标 志 。 但 值得 注意 的 是 , 作者 在 
研究 中 发 现 CO 了 基因 序列 在 按 蚊 种 内 存在 一 定 程 
度 的 差异 , 但 小 于 种 间 差 异 ( 本 组 未 发 表 资 料 ) ， 故 
在 分 析 按 蚊 隐 种 间 的 关系 时 , 需 特别 谨 愤 。 理 论 上 
联合 多 个 分 子 标 志 分 析 可 获得 更 为 符合 实际 的 系统 
发 育 树 ,其 原因 是 信息 量 更 大 。 但 由 于 不 同 的 基因 
具有 相 异 的 进化 速率 ， 故 在 联合 分 析 之 前 , 须 进 行 
数据 的 不 相合 性 检验 ,而 对 于 检验 结果 为 不 相合 性 
时 , 是 否 还 可 进行 联合 分 析 , 各 学 者 持 不 同 观 点 
(Bull et al., 1993; Nixon and Carpenter, 1996; 
DeSalle and Brower, 1997; Hi Ai Al h XR, 2007), 
作者 将 CO IL 5E D3 联合 分 析 时 , ILD 检验 虽 显 示 不 
具 相 合 性 , 但 依据 该 数据 组 重建 的 系统 发 育 树 却 具 
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图 3 依据 “COI + D3” 联 合 数据 组 构建 的 塞 蚊 亚 属 NJ,，MP 与 ML 的 50% 合意 树 
Fig. 3 The 50% consensus tree (NJ, ML and MP) of subgenus Cellia mosquitoes by CO IÍ with D3 conjunction dataset 
分 支 上 数字 为 NJ 树 和 MP 树 的 置信 值 。NJ and MP bootstrap values are above each node. 


有 相当 的 合理 性 。 因 此 , 应 综合 考虑 检验 结 采 、 系 
统 树 的 拓扑 结构 ,以 及 其 他 特征 (如 形态 特征 、 生 
态 习 性 和 传 病 作用 等 ) 判断 联合 分 析 的 可 行 性 。 
依据 形态 特征 研究 按 蚊 属 塞 蚊 亚 属 蚊 种 的 系统 
发 育 关 系 研 究 多 有 报告 (Foley et al., 1998; Sallum 
et al., 2000; Harbach and Kitching，2005 ) ， 主 要 的 
类 群 划分 已 基本 达成 共识 。 而 依据 分 子 特征 的 相关 
研究 目前 针对 系 、 种 团 和 复合 体 成 员 种 的 较 多 ， 
aH: 多 斑 按 蚊 种 团 ( Ma et al., 2006; Walton et al., 
2007) , 众 命 按 蚊 种 团 和 微小 按 蚊 种 团 (Sharpe et 
al., 2000; Garros et al., 2005), 2& & n] Pk iy An. 
sundaicus 复合 体 (Dusfour et al., 2007) LA K 39r SE Iy 
ARF] (Morgan et al., 2009) 等 , 所 及 的 种 类 和 数量 均 
AH BR, 无 法 全 面 甄别 塞 蚊 亚 属 的 系统 发 育 关 系 。 印 
度 学 者 应 用 分 子 特征 构建 了 该 国 塞 蚊 亚 属 部 分 优势 
种 系统 发 育 关 系 ， 以 澄清 形态 近似 种 间 关 系 
(Mohanty et al., 2009) 。 本 研究 涉及 中 国 塞 蚊 亚 属 
蚊 类 有 20 种 , 具有 较 强 的 代表 性 , 有 助 于 评价 现存 
的 按 蚊 分 类 系统 , 为 进一步 全 面 阐明 蚊 科 的 种 间 亲 
RABE SHB AH E E PE BO 
属 系 间 关系 显示 , HUW A AIA ABA, 其 
间 杀 缘 关 系 近 , 各 分 子 特 征 结果 相互 支持 ; 市 热 旷 
系 中 仅 1 种 ,其 位 置 不 稳定 ,正确 与 否 难 以 判断 ; 
新 塞 蚊 系 在 依据 CO 序列 构建 的 系统 树 中 为 并 系 ， 
在 D3 序列 和 联合 数据 组 中 却 为 单 系 , 但 后 者 的 置 


fa EA EK Foley 2# (1998) 和 Morgan 2# ( 2009 ) 也 
支持 新 塞 蚊 系 为 并 系 , 但 各 自 依 据 的 分 子 特征 和 蚁 
种 不 同 , 故 新 塞 蚁 系 是 否 存 在 分 化 需 增 加 种 类 进 一 
步 确 定 。 

蚊 种 的 准确 鉴别 是 建立 系统 发 育 关 系 的 关键 ， 
窄 蚊 亚 属 中 有 些 蚊 种 的 形态 变异 较 大 , 尤其 是 种 团 
和 复合 体 隐 种 。 为 避免 由 分 类 带 来 的 误差 , 本 研究 
每 种 均 应 用 多 个 个 体 , 反复 比较 和 核对 , 并 删 去 了 
一 些 无 法 确定 种 特异 序列 的 种 类 ,如 : AREER TE Be 
An. jeyporiensis, 伪 詹 氏 按 蚊 An. pseudojamesi ( 同 物 
异 名 : Rú WE TZ An. ramsayi) Ml 2G RE FH WC An. 
indefinitus 等 。 重 建 的 中 国 塞 蚊 亚 属 种 间 系 统 发 育 
关系 显示 ,大 劣 按 蚊 复 合体 和 微小 按 蚊 亚 群 的 2 种 
杀 绿 关系 均 为 最 近 , SAPD ESTER 
稳定 。 库 态 按 数 在 印度 为 复合 体 , 存在 5 个 隐 种 ， 
日 在 不 同 地 区 传 症 作 用 存在 差异 ，D3 序列 比 对 结 
果 显 示 中 国 分 布 的 应 为 B/C/E 种 (Singh et al., 
2004) , 但 比 对 了 核糖 体 DNA 第 二 内 转录 间隔 区 序 
列 后 (本 组 未 发 表 资 料 ), 发 现 与 库 态 按 蚊 B 种 的 
碱 其 无 差异 (Manonmani et al., 2007) 。 新 塞 蚊 系 的 
多 斑 按 蚁 种 团 在 中 国 分 布 有 5 种 , 拓扑 结构 均 显 示 
除 伪 威 氏 按 蚊 外 的 其 他 4 种 为 单 系 , 伪 威 氏 按 蚁 与 
之 的 亲缘 关系 均 为 最 远 , 且 不 在 一 个 分 文 ,提示 伪 
威 氏 按 蚊 应 重新 归 类 ， 文 持 早 前 的 观点 (Ma et al., 
2006; Walton et al., 2007) 。 本 研究 涉及 的 环 纹 按 
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蚊 种 团 有 3 成 员 种 未 聚 为 1 支 , 但 依据 联合 数据 组 
显示 3 成 员 种 间 关 系 相 近 , 菲律宾 按 蚊 与 雪 足 按 蚊 
的 亲缘 关系 最 近 , 两 蚊 在 形态 上 也 较 难 区 分 。 有 学 
a e 
., 2007) , D3 序列 比 对 结果 显示 中 国 分 布 的 应 为 
CH LF). J RE doc SPE CHE BT 
WIA, 以 及 斯 氏 按 蚊 、 帕 氏 按 蚊 和 美 彩 按 蚊 在 
PEMA PHRASE ae, 其 原因 可 能 是 各 自 
种 团 内 其 他 成 员 种 较 少 所 致 。 
eo, FERRERA 么 完成 或 正在 
T, 海量 的 基因 序列 信息 不 断 被 发 表 , 对 海量 基 
PARR ERRATA IER ren 
的 学 
有 的 研究 正明 系统 发 基因 组 站 在 外交 soo 99 
化 的 基本 问题 方面 有 巨大 湾 力 (Delsuc et al., 
2005) ,应 用 系统 发 育 基因 组 学 重建 蚊 种 的 系统 发 
育 关系 应 是 未 来 研究 的 重要 方 问 。 
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致谢 ”承蒙 云南 省 寄生 上 哇 病 防治 所 顾 云 安 和 王 不 玉 
及 第 二 军医 大 学 杨 曼 尼 和 李 翔 宇 等 协助 采集 样本 ， 
FER TRA VRE ! 
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